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Abstract 

We propose a simple method by which we can explain the magnetic-field dependence 
of the longitudinal ultrasonic attenuation in a type-II superconductor. It gives a 
curve which is in good agreement with experimental data, in paticular, near the 
lower critical field H c \. We compare it with conventional methods, which is not in 
good agreement with the experimental data near H c \ but near the upper critical 
field H c2 . 
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1 Introduction 



The temperature dependence of the ultrasonic attenuation in a superconductor 
is well described by the BCS theory [1]. According to the BCS theory, the 
longitudinal ultrasonic attenuation coefficient in the superconducting state 
relative to that in the normal state is given by 
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The energy gap parameter A(T) depends on the temperature. 
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However, the magnetic-field dependence has not been explained very well. 
Ikushima et al. performed an experiment about the magnetic-field depen- 
dence of the longitudinal ultrasonic attenuation in the mixed state of pure 
niobium [2]. And they explained it theoretically. The spatial average of the 
energy gap, (A(H,T)), was introduced and assumed to be proportional to 
the root mean square of the order parameter. Then (A(H,T)) must be pro- 
portional to the square root of the magnetization, M, near the upper critical 
field, H c2 . They assumed this relation to be valid throughout the mixed state. 
And they calculated the attenuation coefficient in the mixed state, a s , with 
the measured magnetization. The results of calculation give the curve which 
is not in good agreement with the experimental data except near H c2 . 

In this paper, we propose an original simple method to explain the magnetic- 
field dependence of the longitudinal ultrasonic attenuation, and compare it 
with conventional methods. First, we study the original method, which is to 
estimate the ratio of the space ocuppied with vortices. It gives a curve which 
is in good agreement with the experimental data, in paticular, near H cl . We 
also study the conventional method, which is to calculate the order parameter 
by solving the Ginzburg-Landau (GL) equations. It gives a curve which is not 
in good agreement with experimental data near if cl but near H c2 . In each 
method, we use eq. (1), in which {A(H,T)) is used instead of (A(T)). 



2 To estimate the ratio of the space occupied with vortices 

We take 77 to be the ratio of the space ocuppied with vortices. For example, if 
each vortex is assumed to be a cylinder whose radius is the coherence length 
£, i] can be written as 



Here, n is the number of vortices per unit area on the plane perpendicular to 
the magnetic field. We should notice that not only n but also £ varies with 
the field [3,4]. Let us assume that the energy gap vanishes in each vortex and 
equals A(T) in the other area. Then we can introduce the following equation: 



Next, let us introduce a simple and phenomenological assumption: the increase 
in the magnetic-field energy equals that of the free energy when the field 
increases from H cl to H: 
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(A(H,T)) = (l-r ] )A(T). 
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Uh 2 -H 2 c1 ) = Q s (H)-Q s (H c1 ), 
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Fig. 1. The longitudinal ultrasonic attenuation coefficient in the mixed state relative 
to that in the normal state. It is plotted against the applied magnetic field H at 
4.2K. The lower end of the solid curve corresponds to H c i, and the upper end to 
H C 2. The open squares, experimental data, and the dashed curve, theoretical curve, 
are from ref.[2]. The solid curve is obtained from eqs. (1) and (7). 

and the increase of the free energy should be related to rj. Thus we have 



n = A(H 2 - Hi) 



(5) 



Here, A is a constant. Taking into account that (A(H,T)) vanishes at H 
H c2 , from eqs. (3) and (5) we have 
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Combining eqs. (3), (5) and (6), we obtain 
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The solid curve in Fig. 1 is obtained from eqs. (1) and (7). It is in good 
agreement with the experimental data, particularly near H&. 
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3 To calculate the order parameter by solving the GL equations 



Now, we calculate (I*) 2 ). Here, \l/ is the order parameter: \]> = \l/ /e n ; \l/ is the 
order parameter in the absence of field; / and 7 are the normalized magnitude 
and phase of the order parameter. Hao et al. [5] solved the Ginzburg-Landau 
equations with a trial function, 
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Here, and are two variational parameters; r is distance from a vortex 
axis. And, especially for k ~ 5, the parameters were approximated by the 
following formulas: 
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where B = 2tt/(kA cc u) is the averaged magnetic flux density; A cc \\ is the unit 
cell area of the vortex lattice; £„ D is the value of at B — 0. Using eqs. (8), 
and (9) and (10), the order parameter should be given by 
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Here, the integral is taken over one lattice cell, which is approximated by a 
circle centered at a vortex axis and has the same cell area. By calculating 
eq. (11), we find 
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where i? is the radius of the circle, 
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And, in order to write (I*) 2 ) as a function of H, we assume that B is related 
to H by the following equation [5], 
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Fig. 2. The longitudinal ultrasonic attenuation coefficient in the mixed state rel- 
ative to that in the normal state. The solid curve is obtained from eqs. (1) and 
(9)-(15). The dashed curve is obtained from eps. (1) and (7). The open squares are 
experimental data from ref. [2]. 
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Then, if we use the relation, 
(A(H,T)) _ 



2K,£, v K 1 (f 00 £ v ) 

BK^K^Ufl + ZBK) 1 ' 2 



(14) 



A(T) \ 



W) 

Vl>2 ' 

o 



(15) 



The solid curve in Fig. 2 is obtained from eqs. (1) and (9)-(15). It is not in 
good agreement with the experimental data near H c \ but near H C 2. 



4 Discussion 



In this section, we discuss the validity for each of the approximations which 
are adopted to calculate the three theoretical curves. 

The theoretical curve calculated by Ikushima et al. [2] was proposed to explain 
the experimental data particularly near H c2 . They used the approximation 
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that (A(H, T)) should be proportional to the square root of the magnetization 
M. It is valid near H c2 . However, it is not necessarily valid near H cl . 

In Section 2, we introduced simple and phenomenological assumptions and 
approximations. Equations (4) and (5) should be effective near H c i, where 
the number of the vortices is small. However, they may be less effective near 
H c2 . Because the relation between the free energy and r], the ratio of the space 
occupied with vortices, may be more complex when the number of the vortices 
is large. 

The approximation leading to eqs. (9) and (10) in Section 3 is not necessarily 
effective in explaining (A (if, T)), on which the ultrasonic attenuation depends. 
It was originally adopted to explain the magnetization curves in ref. [5]. The 
magnetization should be strongly dependent on the penetrarion depth A. Also, 
the order parameter \1/ in Section 3 should be dependent on A. On the other 
hand, (A(H,T)) should be dependent on the coherence length £. In eq. (15), 
we assumed the behavior of (A(H,T)) to be equivalent to that of (|^|). In 
other words, we assumed the magnetic-field dependence of A to be equivalent 
to that of £. However, the magnetic-field dependence of A is generaly different 
from that of £ [3]. Therefore the magnetic-field dependence of magnetization 
should be generaly different from that of {A(H,T)). Equations (9) and (10) 
may be effective in explaining (A(H,T)) not through the whole mixed-state 
region but only near H c2 . 



5 Conclusion 



The original method discussed in this paper gives the theoretical curve which is 
in good agreement with experimental data particularly near H cl . On the other 
hand, the theoretical curves given by the two conventional methods are not in 
good agreement with experimental data near H c i but H c2 . Though the original 
method is simple and phenomenological, it may give some information on the 
conformation of the vortices in type-II superconductors particularly near H&. 
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